arXiv:1506.02666v3 [gr-qc] 11 Dec 2015 


Brane induced gravity: Ghosts and naturalness 


Ludwig Eglseer,^’* Florian Niedermann,^’ ^ and Robert Schneider^’ ^ 

^Arnold Sommerfeld Center for Theoretical Physics, 
Ludwig-Maximilians-Universitat, Theresienstrafle 31, 80333 Munich, Germany 
^Excellence Cluster Universe, Boltzmannstrafie 2, 851 f8 Garching, Germany 

(Dated: December 14, 2015) 

Linear stability of brane induced gravity in two codimensions on a static pure tension background 
is investigated. The brane is regularized as a ring of finite circumference in extra space. By explicitly 
calculating the vacuum persistence amplitude of the corresponding quantum theory, we show that 
the parameter space is divided into two regions—one corresponding to a stable Minkowski vacuum 
on the brane and one being plagued by ghost instabilities. This analytical result affirms a recent 
nonlinear, but mainly numerical analysis. 

The main result is that the ghost is absent for a sufficiently large brane tension, in perfect agree¬ 
ment with a value expected from a natural effective field theory point of view. Unfortunately, the 
linearly stable parameter regime is either ruled out phenomenologically or becomes unstable for 
nontrivial cosmologies. We argue that super-critical brane backgrounds constitute the remaining 
window of opportunity. 

In the special case of a tensionless brane, we find that the ghost exists for all phenomenologically 
relevant values of the induced gravity scale. Regarding this case, there are contradicting results in 
the literature, and we are able to fully resolve this controversy by explicitly uncovering the errors 
made in the “no-ghost” analysis. Finally, a Hamiltonian analysis generalizes the ghost result to more 
than two codimensions. 

PACS numbers: 04.50.-h, 04.50.Kd, 95.36.-|-x 


I. INTRODUCTION 

General Relativity (GR) viewed as an effective field 
theory (EFT) gives rise to a famous problem: the ob¬ 
served value of the cosmological constant is numerically 
not stable under quantum corrections—an observation 
which is usually referred to as the cosmological constant 
problem [1]. An auspicious way of addressing this prob¬ 
lem relies on modifying GR by weakening gravity on 
cosmological scales and thereby making it insensitive to 
a cosmological constant (for recent reviews of modified 
gravity see [2, 3], and for the so called degravitation 
mechanism see [4-8]). 

A prominent arena for realizing such an infrared modi¬ 
fication is provided by the braneworld scenario (see [9] for 
a review), according to which our universe is a surface— 
the brane—in some higher dimensional spacetime—the 
bulk. The crucial idea then is to absorb the gravitational 
impact of a cosmological constant into extrinsic curva¬ 
ture, invisible to a brane observer [10-12]. In principle, 
such a surface can be realized by a topological defect like 
a domain wall in five, or a vortex in six dimensions. In 
both cases there are known ways of localizing particles 
inside the defect, see [13] and [14], respectively. These 
constructions generically introduce a length scale rp that 
characterizes the transverse width of the defect. How¬ 
ever, as long as we are interested in low energy questions. 
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which take place on length scales £ much larger than the 
defect size, it is not necessary to resolve its microscopic 
details. In particular, for cosmological applications where 
£ is of the order of the Hubble length, £ tq by at least 
30 orders of magnitude. Thus, in this context the defect 
can in good approximation be treated as infinitely thin.^ 
Then, in accordance with the EFT paradigm, it suffices 
to write down all operators that are compatible with the 
symmetries of the theory. In the case at hand this is 
the higher dimensional diffeomorphism invariance in the 
bulk and its four dimensional counterpart on the brane. 
As long as we do not impose any further symmetries (like 
supersymmetry [8]), we find 

S = 5eh + 5big + Sm ■ (1) 

The first term, 

5eh = J d^A - 2A(^)) , (2) 

describes Einstein-Hilbert gravity in H = 4 -|- n infinite 
spacetime dimensions, where is the H-dimensional 
Ricci scalar, the bulk Planck scale is denoted by Mu and 
the bulk cosmological constant by The second term 

is the induced gravity term on a codimension-n brane: 

5big = Ml J d^x - 2AW) . (3) 


^ In higher codimensions, delta sources generically lead to diver¬ 
gences, and therefore have to be regularized, see Sec. II A. How¬ 
ever, as we will see in Sec. IV, some questions can already be 
addressed at the zero-width level. 
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The induced terms are controlled by the induced gravity 
scale^ M 4 and the brane tension A := The last 

term in ( 1 ) is the action for matter fields localized on the 

brane and coupled to the induced metric ; ultimately, 
it contains all Standard Model fields. In its pure form and 
for infinite volume extra dimensions, this theory is known 
as Brane Induced Gravity (BIG) in general, or the Dvali- 
Gabadadze-Porrati (DGP) model in the special case of 
n = 1 [15, 16]. It is important to note that even if we 
tuned the brane induced parameters M 4 and A to zero at 
the classical level, they would be generated nonetheless 
at the quantum level by particle loops on the brane. In 
other words, irrespective of the underlying microscopic 
theory that gave rise to the surface theory, the BIG terms 
have to be included in order to meet the requirements of 
a natural EFT. 

That being said, we arrive at the following puzzle: 
There are several claims in the literature [17, 18] stat¬ 
ing that for A^^^ = A = 0 and n > 2 the theory ( 1 ) 
has a tachyonic ghost in a weakly coupling regime on 
a Minkowski background, thus raising the question: Is 
the EFT description of gravity induced on higher codi- 
mensional surfaces necessarily plagued by instabilities? 
This question is at the core of the present work. An an¬ 
swer was supposedly given in a recent work [19], which 
claimed that the theory is healthy in that parameter 
regime (A = 0), in contradiction to the former results. 

By performing two independent analyses, we are able 
to disprove this claim of [19] and to confirm the tachyon 
and ghost result from the literature [17, 18] for arbitrary 
codimensions > 2. Moreover, we are able to pin down the 
errors that were made in [19] by explicitly correcting the 
old calculations here in Appendix A and fully reconciling 
them with our new analysis. 

However, this is not the end of the story. In [20] and 
[ 21 ] the parameter range was extended to non-zero but 
sub-critical^ values of the brane tension A (but keeping 
for simplicity =0). In [20] a linear perturbation 
analysis was performed and a particular scalar mode of 
the dynamical spectrum of the theory was shown to lead 
to a strong coupling for phenomenologically interesting 
parameter values, which only vanished in a near-critical 
tension regime. In [21] a completely nonlinear but nu¬ 
merical analysis for homogeneous and isotropic matter on 
the brane was performed, showing the existence of both 
a stable and unstable parameter regime. Unfortunately, 
the stable regime had to be ruled out on phenomenolog¬ 
ical grounds. It should be noted that the two works are 


^ For phenomenological reasons, M 4 should be identified with the 
Planck mass. 

® A pure tension codimension-two brane is known to produce a 
conical geometry characterized by a deficit angle in extra space. 
For the critical value of the tension, Ac = 2ttMq, it becomes 
27r leading to a cylindrical geometry, cf. [22] or [21] for a more 
recent discussion. Even larger (super-critical) values lead to a 
different topology and inflating branes [23] ; they are excluded in 
the present work. 


based on slightly different physical assumptions: While 
the analysis in [ 20 ] allowed for small fluctuations of the 
circumference of the brane described by the radion field, 
the one in [ 21 ] assumed a strictly constant circumference, 
or an infinitely heavy radion equivalently, based on some 
underlying stabilization mechanism. To summarize, the 
results of [ 20 ] pointed towards the near-critical regime as 
a remaining window of opportunity, whereas the cosmo¬ 
logical analysis of [ 21 ] excluded all sub-critical parameter 
space but lacked a completely analytic statement. 

The analysis of this work will establish a comprehen¬ 
sive and analytic picture which nicely confirms the re¬ 
sults of [17, 18, 21] and is fully compatible with [20]. 
To be more specific, by focusing on the codimension-two 
model, we solve the full system of bulk vacuum and brane 
matching equations at the linear level^ around a conical 
background geometry introduced in Sec. II. To that end, 
we regularize the setup by wrapping the brane around a 
circle in extra space and thus promoting the brane to a 
codimension-one object. The gravitational effect of the 
brane tension is then completely absorbed by the conical 
defect, while the geometry on the brane is Riemann flat. 
Thus, from the perspective of a brane observer there is 
a landscape of different vacua corresponding to different 
values of A [20] . (For that reason these solutions are inter¬ 
esting with respect to the cosmological constant problem, 
as the brane tension plays the role of a 4D cosmological 
constant.) Like in the approach of [21], we assume the 
existence of a stabilization mechanism in the microscopic 
theory that forces the compact brane dimension to be 
of constant size. On the level of fluctuations, this cor¬ 
responds to setting the radion field to zero, and can be 
achieved by appropriately dialing the angular pressure 
fluctuations. (In Appendix C we drop that assumption. 
Accordingly, the radion field represents an additional dy¬ 
namical degree of freedom. Moreover, in Appendix D we 
provide an explicit example of a dynamical stabilization 
mechanism.) 

In Sec. Ill we introduce the vacuum persistence am¬ 
plitude (0|0)t as a diagnostic tool, which enables us to 
probe for the presence of a ghost in the fluctuation the¬ 
ory. This is based on the fact that a ghost mode, treated 
with the standard Feynman prescription, would lead to 
a transition probability |(0|0)t|^ larger than one. 

A Hamiltonian analysis in Sec. IV offers a different but 
compatible perspective and generalizes the results to ar¬ 
bitrary codimensions for vanishing A. It also enables us 
to do a solid counting of degrees of freedom in arbitrary 
dimensions. Furthermore, it corrects the erroneous ex¬ 
pression for the Hamiltonian derived in [19]. The final 


^ In general, the linear analysis is sufficient for addressing weak 
field questions like e.g. solar system predictions. Note how¬ 
ever that for certain parameter choices there might emerge a 
Vainshtein-like regime around massive sources [24]. Therefore, 
in this work we limit ourselves to sources that can be described 
within a weakly coupled theory. Moreover, for cosmological ques¬ 
tions one needs to resort to a nonlinear analysis as in [ 21 ]. 
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discussion in Sec. V confirms the consistency of (1) as an 
EFT and thereby answers the initiatory question with 
“no”: the instability in codimension two BIG disappears 
if the brane tension is not tuned unnaturally small. 


Summary 

To summarize, our investigations lead to the following 

main findings: 

• We derive a function /(Mg, M 4 , A, tq) of the model 
parameters, the sign of which determines the healthy 
(stable and ghost-free) and pathological (unstable and 
ghostly) regions in the (sub-critical) parameter space. 
This function is the static version of the one found in 
the cosmology analysis of [ 21 ]. 

• For a natural choice of brane induced parameters the 
theory is in the healthy regime. To be specific, there are 
stable on-brane Minkowski vacua, for which the tension 
is completely absorbed into extrinsic curvature. The 
instability only emerges if the induced gravity scale M 4 
is large while the tension is kept small. 

• The stable vacua are ruled out phenomenologically ex¬ 
cept for the near-critical regime. However, the cosmo¬ 
logical analysis in [21] showed that for FLRW symme¬ 
tries on the brane, i.e. a background geometry with 
non-vanishing 4D Hubble parameter H, the stable 
regime becomes significantly smaller, thereby closing 
the remaining window of opportunity (at least for sub- 
critical brane tensions). 

• The ghost found in [17, 18], the strong coupling ob¬ 
served around static sources in [ 20 ] and the instability 
encountered in [ 21 ] is produced by the same scalar de¬ 
gree of freedom s and thus the manifestation of the 
same pathology: for model parameters with / > 0 the 
scalar mode s is a tachyonic ghost with mass m*. A 
priori, the status of the ghost is not clear as it could 
be an artifact of the unstable conical geometry. In¬ 
deed, the background is destabilized due to the tachy¬ 
onic character of the scalar s, which leads to exponen¬ 
tially growing low momentum modes characterized by 
jpj < m*. However, the nonlinear analysis of [21] did 
not show any indication for the formation of a new 
stable background, instead the solutions exhibited a 
run-away behavior. 

• Higher codimensional versions share the same pathol¬ 
ogy, at least on the parameter subspace characterized 
by A = 0, in accordance with [17, 18]. 

• For more than one codimension, the theory has six 
gravitational degrees of freedom that couple to on- 
brane sources and are thus phenomenologically rele¬ 


vant® . The total number of degrees of freedom is found 
to be D{D — 3)/2 which is the same as in Einstein grav¬ 
ity. 

• Without external stabilization, the radion is tachyonic 
in a broad (and relevant) parameter regime; it thus 
destabilizes the background geometry. In other words, 
assuming some sort of underlying modulus stabilization 
is required for having a consistent regularization. 

Conventions 

As our sign convention we adopt “-I- -|- -|-” as defined 
in [25]. We work in units in which c = h = 1. While we 
work in six dimensions in Secs. H and HI, we generalize 
the spacetime dimensions to D = A + n, where n denotes 
the number of codimensions, in Sec. IV and Appendix A. 
The notational conventions for the index ranges are de¬ 
fined according to: 



0, 1, 2, 3, r, (j) 

(in Secs. II & III) 

a, b 

r, (j> 


0, 1, 2, 3, 5, 6 ..., L> 

(in Sec. IV & 
Appendix A) 

a,b,... 

5, 6, ..., D 


0, 1, 2, 3, f 

p,v,p,a 

0, 1, 2, 3 

i,j, ■ ■ ■ 

1, 2, 3 


The index values r and (j) are used to emphasize the use of 
polar coordinates. For symmetrization we use the con¬ 
vention d(iVj) = {diVj + djVi)/2. Four-momentum is 
denoted by p = {oj,p). 

II. COSMIC STRING BACKGROUND 

In this section, we introduce the ring regularization of 
the codimension-two brane and present the correspond¬ 
ing equations of motion, as well as the static deficit angle 
background that is generated by a pure-tension brane. 

A. Regularization 

Just like for a charged string in electrostatics, the grav¬ 
itational field in general diverges logarithmically at the 
position of a codimension-two defect. This problem can 
be taken care of by introducing a transverse width of the 
brane. From a more fundamental perspective we could 
have done this right from the beginning because there 


® If the regularized brane width is not stabilized, there is an ad¬ 
ditional degree of freedom corresponding to fluctuations of its 
angular size. 
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has to be a physical mechanism which creates the brane 
in the first place. Topological defects—as for example 
a Nielsen-Olesen vortex [26] in two codimensions—would 
be obvious candidates. In general, these solutions predict 
a finite transverse spread of the brane. 

Instead, in favor of technical simplicity, we choose to 
describe the brane as a ring iSi of circumference 27rro in 
extra space. Hence, it becomes a codimension-one object, 
separating an interior from an exterior vacuum region. 
This regularization has already proved to be convenient 
in earlier works [20, 21, 27]. Technically, it amounts to 
the replacement® 

5big Ml [ d®x^/^(7^(®)-2A(®)) , (4) 

JMiXSi ^ ' 

where M| = M|/(27rro), = Xj(^■ kt^MI') and is 

the five dimensional induced metric. Of course, at this 
stage the process of introducing a particular regulariza¬ 
tion is arbitrary. Instead, we could have smeared the field 
over a disc of radius rg or similarly used a blurring func¬ 
tion [18]. We could have even resolved the brane in more 
fundamental degrees of freedom by using the aforemen¬ 
tioned Nielsen-Olesen construction. However, the crucial 
point is that we limit the range of applicability of the 
theory to scales I tq. In other words, as long as we 
are interested in low energy questions, as they arise for 
example in the context of late time cosmology, we do not 
have to resolve the microscopic brane physics; in partic¬ 
ular, we expect low energy questions not to depend on 
the chosen regularization scheme. 

In order to have a consistent regularization, we need 
to stabilize the circumference of the ring. Effectively, 
this can be realized by imposing a certain (non-constant) 
amount of external pressure in angular direction which 
prevents the circumference of the brane from fluctuating. 
In a fundamental picture, this corresponds to a model 
where the angular mode, describing fluctuations of the 
circumference, is very heavy compared to the typical en¬ 
ergy scale . In Appendix D we discuss an explicit ex¬ 
ample of such a dynamical stabilization mechanism and 
show that our effective method is appropriate. 

In that context, let us note that there is a lower bound 
on the mass of the angular mode (or radion)—as well 
as the inverse regularization width 1 /rg—on the basis of 
post Cavendish experiments which is ~ 10“®eV, corre¬ 
sponding to a Compton wavelength of ~ lOO/rm (for a 
recent work, see [29], for a review see [2] and references 
therein). This means that for cosmological applications, 
corresponding to much smaller energies, it is fully justi¬ 
fied to ignore those angular size fluctuations by setting 
the radion to zero or equivalently its mass to infinity. 


® The origin as a codimension-two model requires the existence 
of a regular axis in the interior region, which will be demanded 
throughout this work. Otherwise, the same action could also 
describe a topologically distinct model without an axis [28]. 


However, in order to show that our crucial results are 
independent of that assumption, we also study the case 
where the size modulus is light and thus can be excited 
as an additional degree of freedom, cf. Appendix C. As a 
result, we find that the background geometry gets desta¬ 
bilized as the circumference starts to grow (or decrease) 
exponentially in a broad parameter regime. This shows 
that the stabilization requirement is indeed necessary to 
avoid additional instabilities. The crucial point, however, 
is that the ghost survives in this case, and thus is no ar¬ 
tifact of the ring stabilization. 

By performing a Hamiltonian analysis in Sec. IV (and 
a Lorentz covariant analysis in Appendix Al), we can 
explicitly confirm that the ghost result does not depend 
on the details of the regularization. In that case we do 
not need to solve the equations explicitly and can there¬ 
fore stick to the codimension-two description in terms 
of delta functions. This is an important and successful 
consistency check of our main analysis, which relies on a 
specific regularization. 


B. Equations of motion 

Since the bulk is source-free, the six-dimensional Ein¬ 
stein tensor has to satisfy the vacuum field equations 
away from the brane (r ^ rg), 

G%=Q. (5) 

These have to be supplemented by Israel’s junction con¬ 
ditions [30, 31] across the brane, including the BIG terms, 

- [if“^]) + , (6) 

where 7^“^ is the extrinsic curvature, and square brackets 
denote the jump across the brane, i.e. [/] := /out — /in- 


C. Deficit angle solution 


For a pure tension brane the source reads 





( 7 ) 


where the factor 27rrg is chosen such that A corresponds 
to a 4D brane tension. For the static solution to exist, the 
pressure in angular direction was dialed to zero. In Ap¬ 
pendix D we provide an example how this can be achieved 
microscopically. The exact solution is well known (see 
e.g. [21] and references therein), and can be written in 
the form: 


+ dr^ -j- (/(r)^d(/)^ (8a) 


5(f’) 


r (r < rg) 

^0 + (1 - {r - ro) {r > rg) 


(8b) 
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Here we used a continuous coordinate patch with polar 
coordinates (r, (p) in the two extra dimensions, having 
the standard ranges r e [0, oo), p G [0, 27r). The brane 
is located at r = rg. The geometrical meaning of the 
solution ( 8 ) is that (i) the 4D on-brane geometry is com¬ 
pletely flat and (ii) the 6 D geometry is locally flat as well, 
but has a deficit angle 5 in the exterior, which is related 
to the brane tension by 


Here and henceforth a prime is shorthand notation for dr- 
It is convenient to decompose the 3D spatial components 
of h^B as 

hoi = Ni + d^L , (14a) 

hij = Dij -\- + didjB -\- SijS , (14b) 

hir = -|- diH ', (14c) 


<5 = 


A 

w 


( 9 ) 


In this paper we will only consider sub-critical tensions, 
i.e. A < 2ttMq 5 < 2tt. Otherwise, the static deficit 
angle solution would not be stable, the on-brane geome¬ 
try would instead inflate and the bulk topology would be 
different [23]. 


III. LINEAR STABILITY ANALYSIS 


We will now consider small metric perturbations 
around the deficit angle background ( 8 ), 

gAB = lAB + hAB , (10) 

where 


-fAB = diag (-1,1,1,1,1, g{r)'^) , (11) 

i.e. we still work in polar coordinates 
Note that all indices on first order quantities will be low¬ 
ered and raised with the background metric 'jab and its 

• A R 

inverse 7 . 

The question we want to answer is whether Hab 
contains—at the linear level—instable modes (tachyons 
or ghosts, or both), which can be sourced by an addi¬ 
tional (small) on-brane source i.e. 


-,(5)c 


_ 0/j^(5)c 


1 ^(5)c 


27rro ^ ^ 2 


( 12 a) 

( 12 b) 


Since the source is distributed in a (((-symmetric way, it 
is sufficient to consider metric perturbations that also 
respect this symmetry, because only those are sourced. 
This means that the Hab are (((-independent functions, 
and that = h^r = 0. Furthermore, we can keep the 
brane at the fixed coordinate position rg = const with¬ 
out loss of generality, because a proper motion in radial 
direction (which is not ruled out by stabilizing the p- 
direction), as well as a dependence of the physical brane 
radius on the spatial brane coordinates, can still be ac¬ 
complished by allowing for nonzero components. The 
metric perturbations therefore take the form 


hAB = 




^-N 

hoj 

1' 

0 \ 

hio 

hij 

hir 

0 

1' 

hrj 

hrr 

0 

V 0 

0 

0 



( 13 ) 


where—from a 3D point of view —Dij is a transverse 
traceless tensor, and Ni,Vi, Gi are transverse vectors, i.e. 


D\ = diD^j = 0 , (15a) 

diN^ = diV^ = = 0, (15b) 

while N,l, L, B, S, H are scalars. Even though this de¬ 
composition is not manifestly Lorentz-covariant like the 
approach in [19], it has the great advantage that it is in¬ 
vertible on the space of bounded functions, and so it does 
not introduce any “split ambiguity” (cf. Appendix Al). 
The reason is of course that the Laplace operator Ag has 
no non-trivial bounded solutions, unlike the d’Alembert 
operator ^ 4 . This makes the identification of dynamical 
degrees of freedom much more straightforward. 

The analogous decomposition of the (a&)-components 


hab = VaV&6-|-7abS, (16) 

with Va denoting the covariant derivative with respect 
to the background metric "fab- Explicitly, this gives 

h\ = b" + s, h\=^-b' + s. (17) 

Again, the relation (16) is invertible, because the Laplace 
operator A 2 has an empty kernel. 


A. Gauge-invariant variables 

The linearized bulk theory is invariant under the gauge 
transformations 


5hAB = ^(A^B) ■ (18) 

In order not to spoil the (/>-synimetry, the subject 

to 


e<^ = 0, dA,. = dAr = 0. (19) 

Instead of choosing a particular gauge, we will work with 
a complete set of gauge-invariant variables, which can 


^ Note that in codimension two there is no tensor part, and the 
vector part is absent due to 0-symmetry. 
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chosen to be® 

Aj , s, ( 20 a) 

J:=3S + s, 0:=B + b-2H, (20b) 

P ■-B-b-2{L-I), Q--B-N-2L, (20c) 

A := A - A , W, := 2 A - , ( 20 d) 

where we introduced the dot as shorthand for dt. 

Since we chose coordinates in which the brane is lo¬ 
cated at a fixed coordinate position (r = ro), there is 
a further on-brane restriction on the gauge transforma¬ 
tions, 

C.|o = 0, (21) 

where the subscript “ 0 ” denotes evaluation at the brane. 
This implies that on the brane, there exists an additional 
gauge invariant function, namely 

¥>:=^^lo, (22) 

which will also appear explicitly in the junction condi¬ 
tions below. Physically, it corresponds to the radion held, 
measuring huctuations in the size-modulus of the regu¬ 
larized brane. 


B. Bulk equations of motion 

The bulk vacuum held equations (5) at linear order in 
fiAB read 

c ~ 2VcV(^/i‘"^j 

+ 7 AB = 0, (23) 

where De := These can now be projected onto 

the tensor, vector and scalar components, according to 
the decomposition of the metric perturbations (14). Let 
us emphasize that this projection only requires to di¬ 
vide by Laplace operators, which does not introduce any 
homogeneous functions in the resulting equations. This 
is in contrast to the 4D covariant split adopted in [19], 
where one has to divide by d’Alembert operators, mak¬ 
ing the analysis much more subtle and complicated, see 
Appendix A 1 . 

In the following, we omit all equations which are re¬ 
dundant due to the Bianchi identities. However, since 
we are particularly interested in distinguishing dynam¬ 
ical from constrained degrees of freedom, we only omit 


One might worry that, since the definitions of P^Q and Ci in¬ 
volve time derivatives of some metric functions, one could be 
turning actual dynamical quantities into constrained ones “by 
hand”. However, this is not the case, as can explicitly be seen by 
choosing the gauge B = b = L = Gi = ^ (in the bulk), which is 
always possible. 


those components of (23) which appear in the Bianchi 
identities without time-derivatives. In other words, we 
always keep the stronger equations. Explicitly, we drop 
the and {rr) equations. 

The resulting complete set of bulk equations of motion, 
expressed in terms of the gauge invariant variables (20) 
is: 

• Tensor: 

□eA, = 0 (24) 

This is simply the (zj)-component of (23), projected 
onto the tensor part. It shows that is dynamical, 
carrying two independent degrees of freedom (DOF). 

• Vector: 

AgWi -k 2 O = 0 (25a) 

□eA = 0 (25b) 

The first one is the vector-projected (ir)-bulk equation, 
showing that Wi is constrained. The second one is the 
vector-projection of the (Oj)-bulk equation, with Wi 
eliminated by means of the constraint. Thus, A is 
dynamical, carrying two DOF. 

• Scalar: 

( 2 A 3 -f 3 A 2 ) J -k 4 A 3 S -k 3 A 2 A 3 O = 0 (26a) 

2 j -k A 3 (d -k P) = 0 (26b) 

J - s - Q -k A 3 C) -k P = 0 (26c) 

□ 6 ^ = 0 , 065 = 0 (26d) 

The first equation is the (OO)-component of (23), the 
second one is its (Or)-component (already integrated 
once in r, requiring fall-off conditions in the bulk) and 
the third one is the difference of the {4>4>)- and the 
(rr)-equations (also integrated in r). These three are 
constraint equations that can—for instance—be solved 
for 0,P and Q. Plugging these solutions into the two 
scalar-projected (ij(-components of (23), and taking 
suitable linear combinations, yields the two dynamical 
equations (26d) for J and s. 

In summary, there are 6 dynamical DOF (2 vector, 2 
tensor and 2 scalar), all of which satisfy the 6 D wave 
equation in the bulk. This is the correct number of prop¬ 
agating DOF in six-dimensional GR with azimuthal sym¬ 
metry® . 

Below, we will also see that all of these 6 DOF can in¬ 
deed be sourced. This contradicts the claim in Ref. [19], 
where only 5 sourced DOF were found, allowing to argue 
that the ghost mode (s, see below) would not be dynam¬ 
ical. This wrong conclusion was reached by employing a 
gauge transformation which is in fact not allowed by the 
requirement of SO(2) symmetry, see Appendix A 2. 


® Without this symmetry, there would be 3 additional dynamical 
DOF. 
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C. Junction conditions 


The junction conditions for the two dynamical scalars 
then become 


It remains to derive the linearized junction condi¬ 
tions (6). To this end, it is useful to perform a 3D tensor- 
vector-scalar decomposition of the (perturbation of the) 
energy-momentum tensor, analogous to (14): 

C/o, = , (27a) 

u^3 = . ( 27 b) 


Ml[J'] + M|n4J|o = , 

4M4[s']-2/M|n4s|o = -C/^^, 

where we defined the dimensionless constants 

f -1 _M 

^ ' 4a’ “ ■ 2roM4 ' 


(34a) 

(34b) 


(35) 


Energy conservation then decomposes into 

-Uoo + AsU^^) (28a) 

27rro 

+A3U^^^+u(^’> (28b) 

27rro 

-217P + AsUP = 0, (28c) 

while up and are unconstrained. 

As discussed in more detail in Appendix B, the junc¬ 
tion conditions can be written as 


All the junction conditions, viz. (29), (30) and (34), 
now share the same, DGP-like structure. The only (but 
crucial) difference is that the BIG term in the junction 
condition for s comes with a negative sign if / > 0. 

Note that / matches the function / defined in Ref. [21] 
in the static case (JI = 0), where the sign of / determined 
the stability of the model in the case of FLRW symme¬ 
tries (at the nonlinear level).Therefore, one might al¬ 
ready suspect that the scalar mode s will be a ghost in 
that parameter regime. We will now show on analytical 
grounds that this is indeed the case. 


• Tensor: 

+ MiD4B,,lo = up (29) 

• Vector: 

+ Mia^aio = uP (30) 

• Scalar: 

+ ^5 ^4 J|o = -U'^^ + 3U(^^ , 

(31a) 

4M| (^[V] + + Mi 04 (sjo + 3<p) = 

-U'^^ + 3U‘^^. (31b) 

These last two equations are not yet sufficient to solve 
for J and s, because they contain as a third unknown 
on-brane function. It is determined by the (0(/))-junction 
condition, which can be rewritten as (see Appendix B) 

_ ^ MgroD 4 ((fi - sjo) -k M|D4 s|o = U'*’^ . (32) 

As discussed in Sec. II A, we will now require the proper 
circumference of the regularized brane to be constant, im¬ 
plying ip = Q. (The case without stabilization is discussed 
in Appendix C.) The appropriate U'^^ which is needed to 
achieve this stabilization is determined by Eq. (32), which 
now simplifies to 

U% = {Mi - PMP) □4s|o , := . (33) 


D. Tachyonic ghost 

We will use the vacuum to vacuum transition proba¬ 
bility (in presence of an external source) as a diagnostic 
tool to probe for ghost modes. For the linear theory, it 
is given by 

|(0|0)t|^ = exp[-Im(A)] , (36) 

with 

A-.= r3 (37) 

Here, h“^|o should be evaluated at the classical solution 
in the presence of Ta^. If the probability (36) is larger 
than one (or, equivalently, the imaginary part of A is neg¬ 
ative) then unitarity is violated^^, showing the existence 
of a ghost mode. 


Ref. [21] also uses a second regularization, where the extrinsic 
curvature in the interior is demanded to be static, implying a 
slightly different form of /. We expect—although not explicitly 
studied here—to find the same accordance between both analyses 
in that case. 

As is well known, this unitarity violation does not mean that the 
theory cannot be consistently quantized. Unitarity can indeed be 
restored by choosing a non-standard ie prescription for the ghost 
mode, reversing the sign of the ghost-residue in the propagator. 
However, what cannot be cured is the fact that the Hamiltonian 
is not bounded from below, which—as soon as interactions are 
included—causes a catastrophic instability. This instability is 
already present at the classical level, and has nothing to do with 
quantizing the theory. In any case, a ghost shows that the theory 
is pathological and thus useless. 










Since the ghost mode lies in the scalar sector, we 
can limit ourselves to a a source with vanishing tensor- 
and vector components. Using the bulk equations to 
eliminate all constrained quantities, the source coupling 
term can—using integration by part and energy conser¬ 
vation (28)— then be brought into the form 


_ 7^ 

3 



- 3C7(^)) J|o + U^^s\o 


(38) 


Furthermore, it will be sufficient to consider a source 
satisfying for which only s gets excited 

(i.e. J can consistently be set to zero) and the coupling 
term (37) simply reads 

= '^Jd^x U^^s\o . (39) 

In the following, it will be convenient to work in 4D 
Fourier space, i.e. we introduce 


Im((u) 



Re(G;) 


FIG. 1. Analytic structure of the uj = dependence of the s- 
propagator 1/Z{p), see Eq. (45), for some fixed value |p| 7 ^ 0. 
The poles at ica* correspond to the tachyonic ghost. For 
IpI > m* they lie on the real axis, between the origin and the 
branch cuts starting at the poles at ± |p|; for \p\ < m* they 
lie on the imaginary axis. The contour C (dotted blue line) 
indicates the Feynman-contour of integration. 


s{p,r) := J d'^x e '^'^s{x,r) {p ■ x := pfj,x^). (40) 

The bulk equation (26d) then becomes 

{-p^ -b A 2 ) s = 0 ijP' := Pf,p ^), (41) 

where the covariant 2D Laplace operator with respect 
to the deficit angle background geometry ( 11 ) reads (for 
(()-symmetric fields) 


Ao = di 



dl + -dr 
r 


- /3ro 


(r < ro) 
(r > ro). 


(42) 


The most general solution of (41), which is continuous 
across the brane, regular at the origin, and falls off at 
radial infinity^^, is given by 


s{p,r) 


lojry^) .| 
/o(rov^) 
Aro(Fv^) 
A:o(ro\/?^) ° 


(r < ro) 
(f > ro), 


(43) 


where In and Kn are the modified Bessel functions of the 
first and second kind, respectively, and r := r ^tq. A 
priori, the solution is only defined for > 0. We find its 


Furthermore, for < 0, i.e. for modes which correspond to 
waves propagating in the bulk, one can check that the solu¬ 
tion (43) corresponds to solely outgoing radial waves, if the re¬ 
tarded prescription Im(cj) = -)-e is used, as would be appropriate 
for a classical calculation. This is an important consistency re¬ 
quirement as the brane is the only source of gravitational waves 
in the bulk. Note, however, that below we will use the Feynman 
prescription Im(a;^) = H-e, since we are calculating the vacuum 
amplitude in the quantum theory. 


analytic continuation by choosing the branch cut of the 
square root in the standard way, i.e. along the negative 
real axis. 

Plugging (43) into the junction condition (34b) yields 


4Mi 

-^Zip)s\o = -U^,, 
ro ^ 

with the inverse s-propagator 

^ f 2 f+ I3)z) 
Z[p) ■- afz - z[ —— -b ,, . . 

\Io{z) Ko{{l + /3)z) 

where we introduced the dimensionless variable 

z := ro\/jp 

= rQ\/p^ - . 

The source coupling term (39) finally becomes 


= - 


24M| 


d^p 

(27r)3 




Z{P)’ 


(44) 

(45) 

(46) 

(47) 


where the branch cuts and poles in the w-integration are 
surrounded according to the Feynman prescription, i.e. 
Im (w^) = -b e along the integration contour. The ana¬ 
lytic structure of Z~^, for some fixed value |p| ^ 0 (and 
the case / > 0) is shown in Fig. 1. The branch cuts 
along the real axis can be interpreted as a continuum of 
gap-less Kaluza Klein modes, like in the DGP model [ 6 ]. 
They are also present in pure 6 D GR with a cylindrical 
source, and are thus not expected to cause any problems. 
Below, we will confirm that this is indeed the case. 

However, for / > 0, there are additional isolated poles 
at 


uj = =: iw* , (48) 


where m* is given by m* = z^/ro > 0 , with z* being the 
solution of 


hjz*) Ki{(l + f5)z^) 
io{z*) Kro((i-b/?)z*) ■ 


(49) 
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FIG. 2. Graph of the right-hand side of Eq. (49), which de¬ 
termines the ghost mass m* = z,/rQ. For values of the dehcit 
angle 5 between 0 and 27r, the curve lies in the shaded region. 
It always goes to 2 as a —>■ oo. The gray line corresponds to 
the left-hand side of the equation, for some positive value of 
/. For / < 0 its slope is negative and there is no solution a*. 


The right hand side of this equation is plotted in Fig. 2, 
and shows that there is indeed a solution z* > 0 if and 
only if / > 0. The negative sign of the mass term in 
the dispersion relation (48) shows that this pole in the 
propagator of the scalar mode s is a tachyon, implying 
the existence of exponentially growing solutions for s|o(t). 
Below, we will show that it is also a ghost, in agreement 
with the result in Refs. [17, 18], but generalizing it to a 
background with nonzero deficit angle. 

Even though (49) cannot be solved analytically, one 
can obtain the asymptotic formula for to* in the phys¬ 
ically relevant limit a —)■ oo (i.e. tq <C M|/Mg) by ex¬ 
panding the Bessel function for small arguments, yielding 


2 l-(5/(27r) 

TO*-- (a ^ oo). 

r^jama 


(50) 


Note that for 5 = 0 (and neglecting the small logarith¬ 
mic correction) this agrees with the tachyon mass derived 
in [17, 18], viz. to* ~ Mq/M^. However, the nontriv¬ 
ial deficit angle background leads to an important mod¬ 
ification: as S increases, / approaches zero and Fig. 2 
shows that the tachyon then becomes infinitely heavy, as 
the intersection moves to larger values of z. When the 
threshold / = 0 is crossed, the pole finally disappears 
completely. 

To disentangle the tachyon and branch cut contribu¬ 
tion to we consider two independent integration 

contours in the complex w-plane: a closed path Ci encir¬ 
cling one of the poles, and another open path C 2 running 
along both sides of the branch cut in opposite directions, 
see Fig. 3. It can be shown that the half circle in Ci does 
not contribute to the amplitude^^. It thus follows that 


In Fig. 3 we only show the case when the contour has to be closed 
in the upper half-plane, in which the pole and branch cut on the 
negative real axis contribute. But one can easily check that the 
other case yields exactly the same result. 


/ 





Im(oj) 

A 

V 

\ 


Re(G;) 


^2 


FIG. 3. Decomposition of the Feynman-contour into a closed 
path around one of the poles (Ci) and a branch cut contribu¬ 
tion (C 2 ). The half circle, which closes the contour at infinity, 
does not contribute to 



FIG. 4. Numerical evaluation of the branch-cut contribu¬ 
tion, viz. — Im j^ro J^^dco The positive values im¬ 

ply that there are no ghost modes, irrespective of the sign of 
/. Here, the deficit angle was chosen as 5 = tt, but other 
values do not change this result. 


the original integration along C can be decomposed into 
the sum of the two contributions Ci and C 2 . 

As for the branch cut contour C 2 , we checked numeri¬ 
cally that—at least for an w-independent source—its con¬ 
tribution to the imaginary part of .4^®^ is positive, see 
Fig. 4. Hence, the branch cut contains no ghost modes. 

The contribution of C\ to the amplitude is proportional 
to the sum of the residues of all enclosed poles. Therefore, 
to show that the tachyon—in the parameter region where 
it exists—is also a ghost, let us investigate the residue of 
this pole for the case jpj > to*, i.e. when the poles lie on 
the real axis (as in Fig. 1). (For momenta jpj < to*, the 
pole lies on the imaginary axis and only contributes to the 
real part of which does not affect the vacuum tran¬ 
sition probability (36). Physically speaking, the ghost 
can only be excited for momenta larger than its mass.) 
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A straightforward calculation gives 


-1 

■ ( 51 ) 

Here, the arguments of the Bessel I and K functions 
are z* and (1 + /3)z*, respectively. It turns out that the 
expression in square brackets, when evaluated numeri¬ 
cally, is always positive. However, we did not succeed 
in extracting this information analytically, and therefore 
Fig. 5 shows the contour plot of the residue—leaving out 
the overall factor =Fl/(a;*rQ)—as a function of the two 
independent model parameters <5/(27r) = /3/(l -|- f3) and 
/ = 1 — 3/3/(Aa). For a non-negative, sub-critical tension 
we have S G [0, 27r) and / < 1. Furthermore, as already 
discussed, the ghost pole only exists for / > 0. Thus, the 
plot in Fig. 5 covers the whole relevant parameter space, 
and one can see that the expression in square brackets is 
indeed always positive. 

One can easily check that this result leads to a nega¬ 
tives^ imaginary part of corresponding to a ghost. 
This ghost mode can be excited for 3-momenta |p| larger 
than the ghost mass m*. However, since the ghost is also 
a tachyon, which can be excited with arbitrarily low mo¬ 
menta, the linearized theory is completely unstable (for 
all momenta), if / > 0. 

If, on the other hand, / < 0, the pole and thus the 
tachyonic ghost is absent and the model is stable. The 
stable and unstable regimes are visualized in a parameter 
plot in Fig. 6. It also shows that the tachyon mass di¬ 
verges as the borderline / = 0 is approached. Thus, the 
tachyonic instability is more severe close to the stability 
bound. 


Res 


/ 1 


\zip) 


a; = ± w* I = =F 




2a/ + /3 


-(1+/3) 


E. Regularization independence 

In Sec. HA we argued that we are insensitive to the 
details of the regularization as long as we pose low en¬ 
ergy questions. In the context of the above analysis this 
implies that only ghost poles with a mass m*ro 1 are a 
generic prediction in the sense that they are expected for 
other types of UV models. Contrariwise, this means that 
the precise shape of the stability line in Fig. 6 depends on 
the way of regularizing—with the uncertainty region cor¬ 
responding to the dark shaded red region approximately 
left of the a = 1 line. 

This UV dependence was carefully discussed in [21] 
where two different regularizations were used: 


The sign from the overall factor in (51) is compensated by the 
one from the negative/positive orientation when encircling the 
pole at ittJ*, cf. Fig. 3. Then, there is one more minus sign from 
the explicit overall factor in (47). 



6/2n 


FIG. 5. Contour plot of the ghost residue (51) times 
as a function of the deficit angle 5 and / as defined in (35), 
showing that it is indeed positive in the whole parameter space 
(in which the pole at ai* exists, i.e. for / > 0). The dotted 
lines are lines of constant a, which —> 0 on the left and —> oo 
on the right. 


• One regularization fully resolves the gravitational sec¬ 
tor inside and outside the ring, hence matching exactly 
the regularization used here. A generalized version^® of 
/ valid for FLRW symmetries was derived, which ex¬ 
plicitly depends on the Hubble parameter H. Since H 
is zero for the background solution (8), it would enter 
in (35) only through higher order corrections. The im¬ 
portant point is that by setting H to zero, we exactly 
reproduce the result here. In other words, the old anal¬ 
ysis is in full agreement with the parameter plot shown 
in Fig. 6. 

• Another regularization neglected the interior dynam¬ 
ics. Consequently, the generalized version of / was 
slightly different, in particular, there the borderline 
/ = 0 crossed the (5 = 0 axis, leaving a small stable 
regime near a = 0 in the tensionless case, unlike in 
the other regularization. The corresponding modified 
stability bound is shown as the dashed line in Fig. 6. 
However, it can be checked that the difference is lim¬ 
ited to the regime with m*ro > 1, as expected. In 
particular, this does not affect the phenomenologically 
interesting regime characterized by o: ^ 1, as discussed 
in Sec. VB. 

In summary, we could show that the regularization de¬ 
pendence found in [21] is bound to a regime where the 
inner structure of the brane is already being resolved. 
This had to be expected from an EFT perspective. 


It is denoted by /(t) and discussed in the appendix of Ref. [21]. 
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0.0^—I . I . I..I 

0.001 0.01 0.1 1 10 100 1000 


a 

FIG. 6. Stability of the linearized theory is determined by the 
two model parameters a = Mf/ {2Mgro) and deficit angle S. 
The tachyonic ghost only exists in the red region (/ > 0), 
where the contours indicate its mass (in units of 1/ro). The 
dashed line corresponds to the stability bound in the alterna¬ 
tive regularization of [21]. The dotted line shows the (some¬ 
what stronger) bound on the model parameters discussed be¬ 
low Eq. (70). 

IV. HAMILTONIAN ANALYSIS 

Above, we have shown that a tensionless brane with 
large enough BIG scale is necessarily plagued by a scalar 
ghost mode. In this section, we gather strong evidence 
that the same result holds true in higher dimensions 
{D > 6 ). Thereby, we also provide a complementary 
picture to support the previous result in six dimensions. 
To be concrete, we show that the ghost mode leads to a 
negative contribution to the energy density at the brane 
position. As a diagnostic tool we use the Hamiltonian on 
the constraint surface that we derive for the sourceless 
theory, i.e. Ua /3 = 0. This technique was also used in an 
earlier work [19]. However, there it led to the erroneous 
claim that the theory in 6 D without a brane-tension is 
ghost-free. The reason for that failure is discussed in Ap¬ 
pendix A 2. The purpose here is to present a corrected 
analysis which also extends to higher dimensions. 

We will study small metric perturbations on a bulk 
(and brane) Minkowski background, gAB = Vab + hAB- 
To that end, we use a decomposition of Hab that gener¬ 
alizes the one used in (13), (14c) and (16) to arbitrary 
dimensions. To be precise, we introduce 

hoa =na + dj , (52a) 

hab — dab “f htf^a^b) dadbb Sab^ 5 (52b) 

hia = Eib + diFa + daGi + didaH , (52c) 

where under the SO(n) group dab transforms as a trans¬ 
verse and traceless tensor, En,, Va, Fa, Ua as transverse 


vectors and Gi as a scalar. Moreover, Gi and Fia are 
3D vectors, while all the other fields defined above are 
3D scalars. Note that the definitions for the 3D spatial 
components in (14a) and (14b) still apply. 

Instead of fixing a particular gauge, we will again work 
with gauge invariant variables, i.e. combinations of the 
fields invariant under (18). To that end, we extend the 
definitions in (20) by 

J := 2iS + {n — l)s (53a) 

Fib , dab , (53b) 

Wa •— ‘^Fa Va , Ca •— Via Fa • (53c) 

A crucial benefit of the constraint analysis is that it 
corresponds to a simple reformulation of the theory and 
does not require to find any explicit solution. As a con¬ 
sequence, we do not have to regularize the setup to avoid 
singularities and can instead work with an infinitely thin 
defect described by n-dimensional (5-functions. We will 
check the validity of this approach a posteriori by com¬ 
paring the physical degrees of freedom and their matching 
equations in the special case of n = 2 to the results de¬ 
rived in Sec. HI for a (stabilized) brane when the tension 
is set to zero. We will find perfect agreement. 

Let us stress that this is also a nice way of demonstrat¬ 
ing that our results—in particular statements about the 
dynamical content of the theory—do not depend on our 
way of regularizing the brane as a ring in extra space. 
However, we will see that the two descriptions are only 
equivalent if the (proper) circumference of the ring is de¬ 
manded to be constant. This is clear since without such a 
stabilization the regularization introduces an additional 
degree of freedom corresponding to the angular size mod¬ 
ulus of the brane, which is of course not present for an in¬ 
finitely thin defect. In our notation it is given by the field 
(/?. As discussed in Sec. HA, whether or not we should 
include such a mode depends on the UV model. While 
in the main part of the paper (including this section) it 
is set to zero, we discuss its effects in Appendix C. 

In this section we will use Cartesian coordinates = 
such that the brane is at the constant coordi¬ 
nate position ?/“ = 0. The d’Alembert and Laplace oper¬ 
ator is denoted as before but generalized to D dimensions 
or n codimensions, viz. Dd or A„. 


A. Reduced Lagrangian 

Due to gauge invariance, it is possible to fully express 
the Lagrangian C in terms of gauge invariant variables. 
The calculation is straightforward but little enlighten¬ 
ing; thus, we do not display the rather lengthy resulting 
expression. Varying it with respect to the Lagrange mul¬ 
tipliers P and Q yields the constraint 

2 j -b A 3 (o -f P) = 0, (54) 
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which is identical to (26b) and allows to eliminate P, as 
well as 


M 


D-2 

D 


( 2 A 3 + 3A^)(7 + (2 + n)A 3 s + 3A3A^O 


+ M 2 , 5 „(j/) 2 A 3 J-(n-l)s =0, (55) 


constraining O. 

Before we continue, let us check that the last equa¬ 
tion is indeed a generalization of what we found before. 
Away from the brane, we recover (26a) by setting n = 2. 
Due to the J-term, this equation also contains specific on- 
brane information, which can be compared to the match¬ 
ing equations derived in the last section. This is easily 
achieved by regularizing as described in Sec. II A, which 
amounts to the replacement 

SM - ■ ( 56 ) 

27rro 

Integrating the resulting equation over the interval [rg — 
e, rg-l-e] and identifying 27rroM| = M 4 , one indeed recov¬ 
ers the (00) junction condition^® from Sec. III. This agree¬ 
ment can similarly be checked for all remaining equations 
of motion. 

A third scalar constraint for Q arises in this language 
by differentiating (54) with respect to time and using 
another dynamical equation (which follows from varying 
the action with respect to O) to eliminate O. We find^^ 

3A3Q-2A3J-h3J-(n-I)A3s = 0. (57) 

A similar calculations allows to recover two vector con¬ 
straints, one identical to (25a), and another one con¬ 
straining Wa] for the sake of completeness, 

A 3 W;, + 2C, = 0, (58) 

AnWa + 2Ca = 0. (59) 

We now use the constraints to eliminate the (4 -|- n) 
non-dynamical variables O, P, Q, Wi and Wa (as well as 
their time derivatives) in the Lagrangian, which conse¬ 
quently is expressed solely in terms of dynamical degrees 
of freedom. The resulting Lagrangian is diagonal, and by 
decomposing it into its tensor, vector, and scalar contri¬ 
butions, £ = £t + £v + £sj we get: 


• Tensor: 


4£t = M^-2 - (dAP^jf - (dAdabf 


Mi 


- Sn{y) ■ (60) 


There, the (OO)-junction condition is a certain linear combination 
of Eqs. (34a), (34b) and (B4) with <5 = 0 and </? = 0. 

In the terminology of the Dirac constraint formalism, this would 
correspond to a tertiary constraint (while (54) and (55) are sec¬ 
ondary constraints). 


The 3D tensor describes two DOF which also have ki¬ 
netic support on the brane. They correspond to the 
helicity-two modes in 4D GR and are thus crucial in 
realizing a 4D regime. The extra space tensor car¬ 
ries {n -b I)(n — 2)/2 DOF which all decouple from 
the brane. In particular, they would not couple to a 
localized source on the brane. 

• Vector: 


2£v = M 


D-2 

D 


- {dACif - {dACaf - {dAE^af 


Ml 


-{d^C.f Sniy). (61) 


There are two vector DOF, described by Ci, which have 
a standard DGP-type action and hence a localized ki¬ 
netic term. The remaining 3(n — 1) vector DOF are 
decoupled from the brane. 

• Scalar: 


6£s = 


D-2 

D 


- {dAjy - 


2 («-l)(^ + 2),„ ,2 


■Mi 


2 

^2/ 


-{dAsy 


- {d^jy + {n- lyid^sy S^iy) (62) 


We find two DOF in the scalar sector. They are de¬ 
coupled from each other, and, as expected, s has a 
wrong sign kinetic term on the brane. Also note that 
for n = 1 the scalar s simply disappears, confirming 
that the codimension-one model is indeed ghost-free 
on a Minkowski background. 


In summary, for n > 1 there are always six degrees of 
freedom that can be excited by an on-brane source. Five 
of them {Dij, Xi and J) have a standard DGP-type ac¬ 
tion, whereas the scalar s comes with a wrong sign kinetic 
term on the brane—irrespective of the number of codi¬ 
mensions. Therefore, tuning the tension to zero seems 
to cause the same ghost pathology in any higher dimen¬ 
sion. In addition, there are further [n{n -b 5)/2 — 4] DOF 
that only propagate in the bulk and are invisible to a 
brane observer. In total, we have D{D — 3)/2 DOF cor¬ 
responding exactly to the number of propagating degrees 
of freedom in D-dimensional Einstein gravity. This result 
strongly supports the EFT picture, according to which 
the induced gravity term arises simply by integrating out 
heavy particles in the presence of higher dimensional GR. 
From that perspective it is not surprising that we find 
the same number of propagating degrees of freedom. In 
the codimension-two analysis in Sec. Ill we assumed the 
fluctuation to respect SO(2) symmetry, which effectively 
eliminated three DOF^®. Since the present analysis does 
not rely on that assumption, the result also holds in the 


A similar assumption would here eliminate all degrees of freedom 
that do not couple to on-brane sources. 
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case of additional sources in the bulk, which, in general, 
would spoil the SO(n) symmetry. 

Let us stress that our calculation is also valid for one 
codimension; in that special case s completely drops out 
of the equations, thus leaving a theory of five healthy 
degrees of freedom. 


constitute healthy (brane-coupled) fields. However, the 
expression for the scalar s becomes 


UHs = (n-l)(n + 2)M^-" (H,)" + + (5as) 


- 2{n - If Ml (n,)2 + {d,sf S^{y ), (68) 


B. Reduced Hamiltonian 

A ghost mode necessarily leads to a classical instabil¬ 
ity, which manifests itself in a Hamiltonian that is not 
bounded from below. In the remainder of this section, 
we will show that excitations of the ghost mode indeed 
lower the local energy density at the brane, thereby show¬ 
ing a pathology^®. Since all tensor and vector modes have 
the same DGP-type action as the scalar J, it suffices to 
investigate the scalar sector only. The only difficulty in 
deriving the Hamiltonian consists in a proper treatment 
of the localized terms. This can be consistently done by 
decomposing the conjugate momentum fields 

flj:=dC/6j and ns:=(5£/5s (63) 

into a regular and an irregular part according to 


thus displaying a negative energy contribution at the 
brane position. In accordance with the result of the last 
section, it is present for an arbitrarily small coefficient 
Mi and only vanishes if the induced term is set to zero 
exactly. In the latter case our result simply reflects the 
stability of higher dimensional Einstein gravity. 

As an aside, note that we also derived the Hamilto¬ 
nian from the full Lagrangian and applied the Dirac con¬ 
straint formalism in order to obtain the Hamiltonian on 
the constraint surface. After appropriate redefinitions 
of the conjugate fields, it is possible to check that both 
Hamiltonians are identical. For the sake of simplicity, 
here we only presented the more compact (but equiva¬ 
lent) derivation starting with the reduced Lagrangian. 

V. DISCUSSION 


3nj = 

6fls = 


Mf-^ + MlSn{y)\Uj, (64a) 

{n - l){n + 2) Mf-^ 

- 2(n-1)2 M|(5„(y)ln«, (64b) 


In the first part of this final section, we emphasize that 
our results are compatible with a natural EFT perspec¬ 
tive. In the second part, we discuss the phenomenology 
of the model and offer a brief outlook on future research. 


where we defined 


A. Effective field theory picture 


Hj := j and Hg := s. (65) 

This decomposition is well dehned because J and s are 
regular functions at the brane (only r-derivatives would 
introduce discontinuities). The scalar Hamiltonian den¬ 
sity for both fields, 

Tij := fijj — C,j and T-Ls := LfgS — Cs , (66) 

can then be readily derived. For the healthy scalar we 
get a contribution 


QHj = Mf-^ 


(Hj)" + {dMf + {daJf 


+ Ml 


(Hj)' + {dMf s^{y ), (67) 


which indeed shows that J contributes positively to the 
energy of the system. The tensor- and vector contribu¬ 
tions take the same (manifestly positive) form with J 
replaced by Dij and Ci, respectively. Hence, they also 


Note that we did not calculate the vacuum persistence amplitude 
for n > 2, but we also expect a violation of unitarity. 


The EFT paradigm is based on the assumption that 
low energy physics can be decoupled from high energy 
physics. More precisely, if we are interested in describ¬ 
ing a system at low energies, say below the scale A, we 
do not have to dynamically resolve degrees of freedom 
with masses greater than A. Instead, we write down 
an effective theory in which these heavy particles have 
been integrated out. This is achieved by including in the 
action all operators that respect the symmetries of the 
fundamental theory. In general, those operators come 
equipped with coefficients Ci(A, k) which depend on both 
the cut-off scale as well as the parameters (e.g. masses, 
coupling constants) of the fundamental theory, collec¬ 
tively denoted by k. While the A-dependence is a relict 
of artificially introducing a cut-off and should be taken 
care of by an appropriate renormalization scheme, the 
K-dependence is of physical relevance as it reflects the 
presence of the heavy degrees of freedom in nature. 

When we apply the EFT-reasoning to the BIG model, 
the action (1) should be regarded as the low energy ver¬ 
sion of a more fundamental theory of heavy particles 
that are localized on a four dimensional defect in a D- 
dimensional bulk. By integrating out those heavy par¬ 
ticles, the 4D curvature invariants denoted by iSbig are 
induced on the brane. In this context, the presence of the 
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induced gravity terms is a consequence of the low energy 
description and cannot be avoided in any brane setup. 
Since we know that there are consistent, i.e. ghost-free, 
microscopic theories describing localized particles [14], 
there should also exist ghost-free versions of ( 1 ) within a 
broad regime of induced parameter values. The last qual¬ 
ification is important since we require the parameters of 
the fundamental theory not to be hne-tuned^°. To be 
precise, we expect the process of integrating out heavy 
particles to generate induced operators with coefficients 

A = , Mi = C 2 M* , rg-i = C 3 M* , (69) 

where ci, C 2 and C 3 are dimensionless constants and M* 
is an arbitrary mass scale. We included ro because the 
brane width is ultimately also set by microscopic physics. 
Calculating the exact values for Ci would of course require 
precise knowledge about the fundamental theory. 

For example, in the microscopic model discussed in [14] 
the authors consider a localized scalar particle of mass M 
that couples to the bulk gravity sector. On a perturbative 
level, it can be integrated out by calculating a particle 
loop on the brane. It is now crucial that this loop gives 
rise to both, an induced Einstein Hilbert term as well as 
a brane tension. Moreover, both terms are set by the 
same scale, implying ci ~ C 2 . In particular, setting the 
tension to zero (ci = 0 ) while keeping C 2 nonzero would 
correspond to a fine-tuning, cf. Sec. 7 in [14]. 

Note that while, strictly speaking, the following dis¬ 
cussion only applies to codimension two, we expect the 
qualitative results to be true in any higher codimension 
n > I. When we plug (69) into (35), we obtain 

^ ^ " 2(c 2C3)2 I - J/(27r) ’ 

where the last factor may take any value in the range 
[ 1 , 00 ). 

The main result of Sec. Ill showed that the theory is 
healthy, i.e. ghost-free, for parameters that obey / < 0 . 
Replacing the last factor in (70) by its greatest lower 
bound 1 , we see that a sufficient condition for stability 
is given by (c 2 C 3 )^/cf < 3/2 (shown as a dotted line in 
Fig. 6 ). The important point is that this can always be 
achieved by choosing parameters of order one, i.e. with¬ 
out any fine-tuning. However, once we start to solely 
decrease the parameter ci (in [19] and Section IV it was 
set to zero exactly), the bound gets quickly violated, cf. 
Fig. 6 . In other words, a consistent theory of gravity in- 
dueed on a brane requires the inclusion of a sufficiently 
large brane tension, which is in full accordance with an 
EFT perspective. 


For the sake of simplicity we made the unnatural choice = 0 
for the bulk cosmological constant. However, we will see that it 
already suffices to choose a natural value for the brane tension A 
in order to arrive at a ghost-free fluctuation theory. 


B. Phenomenology 

Independent of the question whether the stable regime 
can be realized without a fine-tuning of parameters, we 
may ask whether it is phenomenologically viable. To this 
end, let us derive the crossover scale of the model in 6 D. It 
can be inferred directly from the brane-to-brane propaga¬ 
tor of the healthy (observable) mode^^ J\o that is being 
convoluted with a constant point source Too oc S^{x). 
The correct expression for the propagator can be ob¬ 
tained from (45) via the replacement / —>■ —2 [to see 
this compare (34a) to (34b)]. The time-independence of 
the source then implies w = 0 and a crossover momen¬ 
tum can be derived by comparing both terms in (45). 
Once they are of the same order, the induced graviton 
dynamics is as important as the bulk dynamics, thereby 
signaling a transition between a 4D and higher dimen¬ 
sional scaling behavior of the gravitational potential. 

To that end, we consider again Fig. 2 which depicts 
the second term of (45) as a function of 2 : = ro|p|. Now 
the crossover momentum can be derived by looking for 
intersections with the linear function 2az. The crossover 
momentum has to be much smaller than the inverse reg¬ 
ularization scale, ro\Pi,\ ^ I; otherwise, we would be sen¬ 
sitive to the unknown UV sector of the theory and—even 
more important—we would enter the higher dimensional 
regime already at microscopic distances below tq which 
is incompatible with observations. By a short inspection 
of Fig. 2, we see that this condition implies a ^ 1 irre¬ 
spective of the value of /3. This means that only the right 
half of Fig. 6 is phenomenologically interesting. While in 
the sub-critical regime this is incompatible with having 
a stable theory, for near-critical values of the tension, i.e. 
1 — S/{2tt) ^ I, there is still a small stable stripe (also 
observed in [ 20 ]) between the stability line described by 
/ and the criticality line d = 2 tt. A lower bound on the 
momentum crossover is found to be (neglecting factors of 
order one) 

ro\Pc\^- fora»l, (71) 

a 

which leads in position space to an upper bound for the 
crossover distance Vf, < arg. This result matches the 
near-critical crossover scale derived in [ 20 ]. 

In order to further assess the phenomenological status 
of the near-critical regime, we have to refer to cosmol¬ 
ogy. To be more specific, instead of considering a static 
(and 4D maximally symmetric background), we have to 
look for solutions with FLRW symmetries on the brane 
which allow for a non-vanishing Hubble parameter H de¬ 
scribing the expansion of the infinite brane dimensions. 
This was done on a fully nonlinear level in [21]. As a re¬ 
sult, a non-static generalization of the function / (again 


Instead, we could also calculate the Newtonian potential for a 
point source which would contain a contribution from s, too. 
However, this would not change the conclusion. 
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separating a stable form an unstable regime) valid for 
FLRW symmetries was derived. The important point is 
that for H ^ 0 the stable near-critical parameter window 
becomes smaller. To be precise, it was found that for a 
stable, sub-critical solution to exist the following bound 
has to be fulfilled^^: 

Hr,<l. (72) 

This shows that the model does not allow for a stable 
4D regime in cosmology (which would demand Hrc ^ 
1). This agrees with the observation in [21] that the 
stable solutions are always governed by 6D dynamics. We 
therefore conclude that the sub-critical and near-critical 
model has to be dismissed in accordance with the analysis 
in [21]. 

Let us stress that the present work and the cosmology 
analysis in [21] complement each other very nicely. In 
both cases we find a stability line in parameter space, 
characterized by a function /, which separates a stable 
from an unstable regime. While in the cosmology analy¬ 
sis the qualification of a certain parameter regime as sta¬ 
ble or unstable could only be inferred by solving the full 
differential system numerically, the present work permits 
an analytic assessment which shows perfect agreement for 
the case H = 0. On the other hand, only the cosmology 
analysis was able to rule out the near-critical parameter 
regime which still seemed to be viable for H — 0. 

Finally, let us comment on the range of applicabil¬ 
ity of the linear analysis presented here. Around static 
and spherically symmetric sources the codimension-one 
model is known to possess a Vainshtein-like radius below 
which the linear approximation breaks down. As it is 
parametrically large compared to the Schwarzschild ra¬ 
dius, it turns out ot be crucial in restoring a 4D regime 
on solar system scales [32-34]. It is thus natural to ask 
whether a similar effect could exist in two codimensions 
as well. Answering this question would require to derive 
higher order corrections to the propagator in (45), which 
is beyond the scope of the present work. Instead, we al¬ 
ways assume to be in a weak coupling regime (which can 
be understood as a requirement on the source terms). 
However, since the (sub-critical) model is already ruled 
out by to the nonlinear analysis in [21], which is obviously 
insensitive to that problem, this issue seems not relevant 
anymore. 

C. Outlook 

So far, it is not clear whether the picture might be more 
favorable in higher dimensions {D > 6). In Sec. IV we 


The formula corresponds to Eq. (50) in [21], Note that there a 
different (/3-independent) definition for the crossover scale was 
used: = 6arQ. In order to allow for the most conservative 

assessment, we express the same inequality in terms of the upper 
bound on the crossover derived here. 


showed that the instability is also present for arbitrary 

> 6 on the parameter subspace A = 0. But at the mo¬ 
ment we do not know how the stability bound depends on 
the brane tension in more than two codimensions. This 
question is left for future research. 

The remaining hope for the model in 6D resides on 
the super-critical regime. In a recent work (without BIG 
terms) it was shown that in this regime there is no sta¬ 
ble Minkowski vacuum on the brane. Instead, the brane 
starts to expand in axial direction at a constant rate [23] . 
With respect to the cosmological constant problem, these 
solutions seem not to help since effectively the tension 
gets shifted by a constant amount (~ M|) which is by 
far to small to help with the problem. Despite their in¬ 
capability of addressing the naturalness question of the 
cosmological constant, they could nevertheless provide 
a viable infrared modification of gravity. We deem the 
super-critical regime therefore an exciting and promising 
playground for future research. 

Another interesting questions concerns braneworld 
models with compact dimensions. This could be stud¬ 
ied within the theory (1) by simply imposing different 
boundary conditions. Consequently, the continuum of 
branch-cut states would be replaced by a discrete tower 
of massive Kaluza-Klein modes, see [35] . Since the ghost 
pole is independent from the cut, we suspect that the 
stability of the model would also be threatened in that 
context. However, an explicit calculation remains neces¬ 
sary. 

Finally, let us recall that so far the bulk cosmological 
constant was always set to zero. It would be inter¬ 
esting to see how the healthy region in parameter space 
gets modified if also a natural value for X^^^ of order M|, 
was included. 
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Appendix A: Correcting the errors of the old 
“no-ghost” analysis 

In the present work we have extended the ghost analy¬ 
sis of codimension-two BIG to a nontrivial brane-tension 
background. As discussed in the main text, the major 
result is that the ghost disappears if the tension is large 
enough, thus reconciling healthy codimension-two BIG 
with a natural EFT expectation. 
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For vanishing brane tension, however, the ghost always 
exists (if the induced gravity scale is nonzero). This result 
is in contradiction to those of Ref. [19], where higher codi- 
mensional BIG was analyzed at the linear level around a 
6 D Minkowski background, and claimed to be stable. In 
this appendix, we resolve this tension by explicitly iden¬ 
tifying the errors made therein. 


1. Lorentz covariant analysis 

One claim in [19] was that the scalar S (see below)— 
which corresponds to the ghost mode identified in earlier 
works (as well as in the present one)—does not threaten 
the stability of the model because it would be a con¬ 
strained quantity. We will now show that this is wrong: 
S is in fact not constrained, but dynamical. 

When we linearize the theory on a 6D Minkowski back¬ 
ground, i.e. for a vanishing deficit angle, the full dynam¬ 
ical system (5) and (6) can be written in a concise way 
by introducing d-functions. In Cartesian coordinates, we 
find 

Sn{y) v) • (Al) 

where is the first order Einstein operator in D 

dimensions 

OmV = + V^^dMdN - 2 5% 

+ r]MN {d^d^ - Dn) , (A2) 

and the 4D localized source is related to the first or¬ 
der energy momentum tensor via tj^i, := —2 We 

decompose the graviton field in a 4D Lorentz covariant 
way, 

-k B + , (A3a) 

hfLb — -\- + d^dbH . (A3b) 

The purely extradimensional components hab are decom¬ 
posed as in (52b). Here, is a transverse and traceless 
4D tensor, V^, £^t,, Gfj, are transverse 4D vectors and B, 
S, J-b, H as well as all functions appearing in (52b) are 
4D scalars. Furthermore, J-b and 5^6 transform as vectors 
under the SO(n) group. 

Even though this decomposition makes the Lorentz co- 
variance of the dynamical equations manifest, it has the 
general disadvantage that the components are not deter¬ 
mined uniquely. More precisely, this ambiguity can be 
parametrized in terms of a set of homogeneous functions 
where here and henceforth “homogeneous” refers 
to solutions of the 4D homogeneous wave equation, i.e. 

= 0. It can be easily checked that the decompo¬ 
sition (A3) is then invariant under the transformations 


II 

"<3 

(A4a) 

SB = x^^^ - Ax^^^ , 

U 4 

(A4b) 

sv, = x!?> + ^a„x“>, 

(A4c) 

II 

4 9^5^, A _ 3 3^^ 

xW 



(A4d) 

1 as 



sn = x^"^ , STa = xi"^ , 

(A5a) 

5Gb. = - 


(A5b) 


where and are subject to the two conditions 

= -35^ , a“xi®^ = 0. (A6) 

Here and henceforth, (l/n 4)'0 is a shorthand notation for 
the convolution of tp with the retarded Green’s function 
of the 4D d’Alembert operator ^ 4 . 

Besides the split ambiguity, there is the usual gauge 
freedom (18). Instead of fixing a particular gauge, we 
will again use a complete set of gauge invariant variables, 
viz. S, s,'Df^^,£fj^b and dab, as well as 

0:=B + b-2n, (A7a) 

Wfa := Gf.-V^, (A7b) 

ya ■■= Ta-Va- (A7c) 

Although these are invariant under gauge transforma¬ 
tions, they are not invariant under the homogeneous 
transformations (A4) and (A5), but transform as 

dO = -^X^'^+X^'^-2x^"^ (A 8 a) 

U 4 

- x|f^ , (A 8 b) 

= xi'^ • (A8 c) 

Now it is crucial to realize that one certain combination of 
the functions x^^^) xl^^ and x^'^^ does not affect the gauge 
invariant quantities. Explicitly, one can easily check that 
they only enter via the combinations^^ 

:= *<“> - 2x<-‘>, if := xf + a„x'“ , (A9) 

(or any linear combination thereof) in the relevant 
Eqs. (A4d), (A 8 a) and (A 8 b). 


Note that is subject to the same relation (A6) as ^rid 

thus also only has three independent components. 
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In summary, the gauge invariant variables 
{S,s,0,yV^,yaj'D^^,S^b,dab} are only unique up 
to the 4 + n independent homogeneous functions 
Xa^^}- This will be crucial for correctly 
inferring the number of dynamical degrees of freedom 
(DOF), because each of the x’s can be used to eliminate 
one would-be dynamical component. 

To derive the equations of motion for our variables, we 
start by investigating the /r6-components of (Al). Taking 
their double divergence d^db leads to 

□ 4 A„ (35-b (n-l)s) = 0 . (AlO) 

Demanding fall-off conditions at spatial infinity allows 
to simply drop the extra space Laplace operator A„. 
The general solution then becomes 3iS-|-(u—l)s = 
with an arbitrary homogeneous function. Instead of 
choosing initial conditions to fix we make use of the 
split ambiguity (A4a) parametrized by to remove 
it from the equation, s is therefore constrained by the 
relation 

s =-^5 , (All) 

n — 1 

and therefore no independent DOF. 

By acting with a single divergence or db on the /r6- 
components of (Al) and using (All), we obtain 

O^AnWfj, = 0, □4A„3^tj = 0, (A12) 

respectively. As before, A„ can be dropped, and the 

f 3) 

freedom to choose initial conditions gets “eaten” by x)i 
and Xa^\ yielding 

= 3^a = 0. (A13) 

These relations can be used to simplify the fib-sector of 
(Al), leading to the wave equation 

DoSt^b = 0. (A14) 

Thus, the fields £^b constitute 3(n — 1) DOF that are not 
sourced by brane-localized matter fields. 

There is only one freedom left in choosing the decom¬ 
position, namely the function We will use it in the 
same way as before to derive a constraint equation for O. 
After taking the trace of the a&-components of (Al) and 
using (All), we find 


Once we plug this solution back into the a&-components 
of (Al), we find a wave equation for the transverse and 
traceless SO(n)-tensor modes, 

nDdab = 0, (A17) 

which therefore constitute further [(n-|-l)(n —1)/2] DOF 
that are not coupled to on-brane matter. 

Finally, let us consider the /xi/-components of (Al). 
Taking its trace and making use of all solutions we de¬ 
rived before yields an equation for the scalar 5, 

^ + 6M| 045 ) <5„ {y ). (A18) 

Since we already made use of all shift ambiguities, S is 
a real dynamical mode. In fact, for n = 2 it coincides with 
the tachyonic ghost mode from the analysis in Sec. Ill 
and hence fulfills the same dynamical equation. In order 
to demonstrate the equivalence, we use (All) to rewrite 
(A18) in terms of s and use the regularization introduced 
in Sec. IIA by simply replacing 

hiy) - ‘AplA . (A 19 ) 

27rro 

Integrating the resulting equation over the interval [ro — 
e, ro -I- e] yields the matching equation (34b) with 5 = 0 
and = C/'^^/(27rro); the vacuum equation (26d) is 
obtained by evaluating the equation for r ^ tq. 

Having thus established the ghost character of S (or s 
equivalently) for n = 2, we turn to the remaining sourced 
DOF in the general case. In order to derive the dynamical 
equation for the 4D tensor we need to solve for S 
more explicitly. To be precise, S fulfills the equation 

where we introduced the homogeneous function 
which keeps track of the freedom to choose initial condi¬ 
tions for S. (Recall that there is no split ambiguity left 
which could be used to set to zero. This subtlety 
was missed in [19], cf. the discussion below.) This allows 
to derive an equation for the tensor 


□4O = - 


n — 2 
n — 1 


S . 


(A15) 


A priori, this is a dynamical equation for O sourced by 
S. However, according to (A8a) the decomposition is 
invariant under the shift SO = x^^^ • Again, this implies 
that O is not a true DOF as we can impose arbitrary 
initial conditions without affecting the solution for /iab- 
For instance, we can choose x*'^^ such that O becomes 


o = - 


n-2 1 ^ 
n — 1 ^4 


(A16) 


- Ml □ 4 I?;..) 5„(j/) 

-I- , (A21) 

with the transverse and traceless tensor U^i,' defined as 

^ • (A22) 

This equation is not completely decoupled, since it de¬ 
pends on the initial conditions for S through Ac¬ 

cording to the analysis in Sec. HI, the sign of the induced 
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term does not imply a ghost. As is also subject to the 
five constraint equations = 0 and = 0, it de¬ 

scribes five healthy DOF. Therefore, the theory contains 
6 sourced {5, and [(3-|-n/2)(n— 1) — 1] non-sourced 
{Sfj,b, dab} DOF, which makes a total of D{D — 3)/2, cor¬ 
responding to the number of DOF in D-dimensional OR, 
in accordance with the results of the Hamiltonian analy¬ 
sis in Sec. IV. 

To summarize, the manifestly covariant analysis ex¬ 
actly recovers the results of the analysis in Sec. Ill and 
generalizes it to arbitrary higher codimensions. For 
n — 2, the five healthy DOF can be identified accord¬ 
ing to —>• {Dij,Ci,J}, whereas in both cases the 

ghost is described by S (or s equivalently). 

In the remainder of this appendix, we want to point out 
the error in the analysis of Ref. [19]. There, the identi¬ 
cal dynamical equation (A18) for S was found, indicating 
that iS is a ghost. It was then argued that there would 
be a further, more restrictive constraint equation for S 
which could be derived from the (OO)-component of the 
modified Einstein equations, rendering S non-dynamical. 
For its derivation it was necessary to use of the dynami¬ 
cal equation for However, that equation in Ref. [19] 
differs from (A21) by the last term The rea¬ 

son is that the equation was derived by applying the 
transverse-traceless projector (for a definition 

see the appendix of [19]) to the equations of motion (Al) 
and tacitly assuming that it commutes with the linearized 
Einstein operator ■ However, in general the com¬ 

mutator of those operators is non-zero but yields a ho¬ 
mogeneous function due to the occurrence of the Green’s 
function 1/04 in the projector. This fact caused the fail¬ 
ure of the analysis in [19]. 

More explicitly, for an arbitrary function '0 the follow¬ 
ing relation can be derived easily: 


where are four homogeneous functions that are 
uniquely determined by ^/>. This expression could be 
generalized to the commutator of the projector 
with 0^1^^ Instead of applying the projection oper¬ 
ator explicitly, we decided to derive (A2I) by consecu¬ 
tively solving the scalar and vector sector and keeping 
track of all homogeneous functions as described above. 
By doing so, we avoid all potential pitfalls related to 
the non-commutativity of the projector. Finally, it is 
straightforward to check that the (OO)-Einstein equation 
is identically fulfilled and does not provide any new in¬ 
formations; in particular, it is no constraint on S. More 
explicitly, we can eliminate all fields except for S from it 
and arrive at an expression which is equivalent to equa¬ 
tion (A20). However, due to the occurrence of the ho¬ 
mogeneous function (which was missed in the old 
analysis), this should not be interpreted as a constraint 
on iS, it is simply a different (but equivalent) formulation 
of the dynamical equation (AI8). 


= (A23) 


2. Hamiltonian analysis 

The analysis of Sec. IV showed that the Hamiltonian 
on a Minkowski background for n > I is not bounded 
from below, which is another manifestation of the ghost 
mode s. This is another clear contradiction to the re¬ 
sult in [19], where it was claimed that the Hamiltonian 
is positive definite for n = 2. To resolve it, we have 
to reconsider the gauge choices in combination with the 
symmetry assumptions that were made in [19]. 

Beside the choice H- = 0, the following two gauges 
were imposed, cf. Eqs. (47)-(49)^^ in [19], 

dihij = 0, /i56 = 0. (A24) 

In order to check whether these gauge choices are con¬ 
sistent and leave any residual gauge freedom, we have to 
infer their transformation behavior according to (18), 

i (Aa^, + , (A25) 

= 9(5^6) ■ (A26) 

It follows that the first gauge completely fixes without 
any residual gauge freedom. The second gauge leaves 
a residual gauge freedom in the ^a“Sector which is sub¬ 
ject to the condition 9(5^6) = 0- Furthermore, the gauge 
functions have to respect the SO(2) symmetry that 
was also assumed in [19]. Explicitly, this requires = 0 
and d(j,^r = 0, which in Cartesian coordinates implies 
^[ 5 ^ 6 ] = 0- This relation together with the residual 
gauge condition completely fixes the two functions 
and does not allow for any further gauges. In partic¬ 
ular, it is no longer possible to gauge the combination 
2drdihGi — cos(()))(A 3 -I- A 2 )/i 55 as it was done in Eq. (59) 
in [19]. However, this gauge was crucial in that case 
to demonstrate the positive definiteness of the Hamilto¬ 
nian. Correspondingly, the corrected constraint analysis 
of Sec. IV comes to a different result. 


Appendix B: Junction conditions 


Here, we present the derivation of the junction condi¬ 
tions (29)-(32) in more detail. A straight-forward calcu¬ 
lation yields the following non-vanishing components of 
the extrinsic curvature tensor at linear order around 
the deficit angle background (11): 


^ {drh\ - d^Kr - d,h\)\^ , 


(Bla) 

(Bib) 


All equation numbers here refer to the Journal version of 
Ref. [19]. 
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The junction conditions are all obtained by plugging this 
(and the linearized 5D Einstein tensor) into ( 6 ), pro¬ 
jecting onto the desired components, and simplifying 
the BIG terms by means of the bulk equations in the 
limit r —>■ rg. For the tensor and vector 

components, this immediately gives the junction condi¬ 
tions (29) and (30). Similarly, the two scalar junction 
conditions (31) are obtained by taking the required lin¬ 
ear combinations of the 4D trace, scalar and 

equations. 

Next, we have to consider the (0^) component. Af¬ 
ter using the jump of the r-derivative of the bulk equa¬ 
tion (26c), it takes the simple form 


The two junction conditions (31b) and (32) can still be 
used to derive a closed equation^® for s. 


Mils': 


■Ml 



□4 


P-2a 


2a (1 -I- P) rl_ 


s|o 



(C2) 


After performing a 4D Fourier transform and using the 
general bulk solution (43) for s, we arrive at 


4Mi ~ 

-MLz{p)s\, = 
ro ^ 


(C3) 


Mla^lb'] + Mla^slo = u \. (B2) 

Furthermore, the second equation in (17) shows that con¬ 
tinuity of implies 


where now the inverse s-propagator is given by (in terms 
of 2 : := ro\/p 2 ) 


Z{p) := afiz^ + f2-z( 


Kiiil+(3)z) \ 

Ko{{l+l3)z)J 


, (C4) 


[b'] 


s|o) 


rpS 

211 — 5 


s|o) ■ 


(B3) 


Using (B3) in (B2) yields Eq. (32). 

For completeness, let us note that the jump in the r- 
derivative of the constrained scalar O is determined by 
the component of Israel’s junction conditions. Af¬ 

ter using again the bulk equations, as well as (B3) and 
the fact that [H'] = 0 (due to continuity of the metric), 
it simplifies to 


Mi\[0'] + 


rpS 

271 — 5 


Ml 


(s|o - ‘f) 


J — s 


— Q ^ 


= , (B4) 


with 


h :=2 



_ 13-2a 

1 + /3 


(C5) 


This is very similar to the inverse s-propagator in the 
stabilized case, Eq. (45). The only difference is that the 
coefhcient / is slightly modified into fi and—more im¬ 
portantly, as we will see below—there is an additional 
constant (i.e. p-independent) term / 2 . 

As before, we can restrict ourselves to sources for which 
only the scalar mode s is excited. This can be achieved 
by setting all tensor- and vector source terms to zero, and 
requiring 





(C 6 ) 


where an overall A 3 was dropped. The jumps in 
the r-derivative of all the remaining constrained gauge- 
invariant quantities (Wi,P and Q) can readily be ob¬ 
tained from the r-derivatives of the corresponding bulk 
equations. 

As a consistency check, we explicitly verified that all 
the junction conditions, together with the bulk equations, 
imply the energy conservation equations (28), as is guar¬ 
anteed by the Gauss-Codazzi equations. 


In this case, the scalar mode J can again be set to zero, 
and the full source vertex (37) takes the same form as 
before (47), with the replacement Z ^ Z. Therefore, 
the stability analysis is completely analogous. 

It can now easily be checked that the propagator Z~"^ 
has a tachyonic pole (<t4> Z is zero for some 2 * > 0) if and 
only if /i > 0 or /2 > 0. The corresponding regions in 
parameter space are depicted in Fig. 7. They are disjoint, 
and separated by a narrow (but finite) stripe in which the 
tachyon is absent and the model is thus linearly stable. 


Appendix C: Non-stabilized circumference 

In this appendix, we investigate the case when the 
brane circumference is not kept constant. This simply 
means that the radion p is not set to zero, and is 
arbitrary. For simplicity, we will only consider the case 

(Cl) 


Note that when comparing this equation to the one derived in 
Ref. [20], viz. Eq. (5.41) therein (the corresponding scalar mode 
is called X in [20], and is related to ours via s = —3X), one has to 
take into account that the energy momentum tensor in [20] 
differs from our is defined in the Einstein frame and 

thus satisfies standard energy conservation {d^r^jy = 0), whereas 
is defined in the Jordan frame, for which, on the deficit angle 
background, additional terms cc A have to be included, cf. our 
Eq. (28). We thank Nemanja Kaloper for clarifying this point. 


ut = o. 
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The case /i > 0 is very similar to the condition / > 0 
in the stabilized case. Indeed, the delimiting line / = 0 
in parameter space, shown as a dashed line in Fig. 7, only 
gets shifted by a small amount. 

The case /2 > 0, however, implies that there is also a 
tachyon in the upper left region in parameter space. At 
first, it might be surprising that this region is also unsta¬ 
ble, because it also includes the case a = 0, which is just 
pure 6D GR without any induced terms and should be 
a healthy theory. The resolution to this puzzle is rather 
simple: Evaluating the residue of Z~^ at the tachyon pole 
as we did in Sec. Ill D, we find that, while it is again neg¬ 
ative for fi > 0, it is positive for /2 > 0. In other words, 
the tachyon is only a ghost in the lower right region. In 
the upper left region, the tachyon is not a pathology, but 
merely a reflection of the fact that, without fixing the 
brane circumference, the static deficit angle background 
is not stable. Instead, the brane wants to expand (or 
collapse) in radial direction. 

To summarize, the ghost-criterion is basically indepen¬ 
dent of whether the circumference is stabilized or not. In 
particular, the naturalness discussion from Sec. V A still 
applies to the case of free radial expansion. On the other 
hand, the static deficit angle solution is not stable un¬ 
der angular size fluctuations. To overcome this problem, 
we have to make additional assumptions about the un¬ 
derlying microscopic model. In fact, from a fundamental 
perspective, the existence of some sort of stabilization 
mechanism has to be expected as there are known sta¬ 
ble vortex configurations in two codimensions, cf. [26]. 
Fixing the proper circumference turned out to be a con¬ 
venient way of realizing such a mechanism in an effective 
low energy description. In the following appendix we 
study an explicit microscopic example. 


Appendix D: Example of a stabilization mechanism 

In this appendix, we will present an explicit realization 
of a stabilization mechanism, which is capable of keeping 
the circumference of the ring-regularized brane (approxi¬ 
mately) constant. In the main text, it was assumed that 
some mechanism of this kind could exist, and so it is in¬ 
structive to see a concrete example. Following [20], we 
put a massless scalar field F on the brane, i.e. we add 

= (Dl) 

to the action (4). We then choose a background solution 
for which the scalar field winds around the ring, 

E = g<^, (D2) 



a 


FIG. 7. Without stabilization (and = 0), there are three 
regions in parameter space: in the lower right region, where 
/i > 0, there is a tachyonic ghost; the delimiting line is almost 
the same as in the non-stabilized case, shown as a dashed line. 
In the dark green region in the middle /i ,/2 < 0, and the 
model is linearly stable. In the upper left region /2 > 0 and 
there is a tachyon which is not a ghost. It shows that the 
static deficit angle background is not stable if is not used 
to fix the brane circumference. 

where g is a constant. Together with (4), this leads to the 
following background energy momentum tensor (EMT), 

-b Si !;, (D3a) 

, (D3b) 

where gSS = l/rp at the background level. As discussed 
in Sec. IIC, the static pure tension solution exists for 
= 0, which can now be achieved by choosing q 
such that 

^ . (D4) 

4r2 

From (D3a) we see that this doubles the contribution of 
A*-^^ to the 4D EMT. After absorbing this factor by a 
trivial renormalization^®, and identify¬ 

ing the 4D brane tension via 

A = 27rroM|A;®]„) , (D5) 

we arrive at the background EMT (7). 


In the main text, we work with the renormalized quantity 
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So winding E around the ring in this way indeed leads 
to the static deficit angle solution presented in Sec. IIC. 
But a successful stabilization should also suppress the 
fluctuations of the brane circumference, as measured by 
the radion field (p = Let us now show that this can 
also be achieved in this particular example. 

First, note that the fluctuations of E can be consis¬ 
tently set to zero, since this field is not sourced. Next, 
perturbing the background metric in (D3) leads to the 
following first order contribution to the EMT, 



Mj S 

To 27r 




(D6) 


where we eliminated in favor of the deficit angel 
6 . [There are similar contributions to which cancel 
the corresponding (^-terms on the left hand side of (31).] 
Plugging this into (32), it becomes a mass term for the 
radion. 


(□4 - ml) (p + 



, 1 - 5/2tt 

^- 


' 0 


(D7) 


Since oc l/rp, the radion gets heavier as the regular¬ 
ization size ro decreases, and so it costs more energy to 
excite this DOF. In the low energy regime we are aiming 
at, i.e. at energies well below 1/ro, the kinetic term of p 
is negligible,and we obtain 

/ 2 a \ 04 

The remaining explicit p-tevva in (31b) is thus negligible 
compared to sjo and can be dropped, while can be 
replaced by 


U%^{Ml-pM^ro)U4s\^. (D9) 

This is exactly the same as (33), and we indeed recover 
the scalar equations (34) from the main text. While 
there, they were derived by assuming some underlying 
stabilization mechanism, and then inferring the required 
from the field equations, we have now seen an explicit 
example that gives the radion a large mass, and consis¬ 
tently yields the same and field equations after the 
radion has been integrated out. 
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